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We propose a definition for topological order at nonzero temperature in analogy to the usual zero 
temperature definition that a state is topologically ordered, or "nontrivial" , if it cannot be trans- 
formed into a product state (or a state close to a product state) using a local (or approximately local) 
quantum circuit. We prove that any two dimensional Hamiltonian which is a sum of commuting 
local terms is not topologically ordered at T > 0. We show that such trivial states cannot be used to 
store quantum information using certain stringlike operators. This definition is not too restrictive, 
however, as the four dimensional toric code does have a nontrivial phase at nonzero temperature. 
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Topological quantum computation^ is one of the most 
promising approaches for building a quantum computer. 
The nonlocal encoding of quantum information in the 
ground state subspace of topologically ordered systems 
protects it against decoherence. Topological order is also 
one of the most interesting current topics in condensed 
matter physics, as models such as the Levin- Wen models^ 
provide a method for constructing novel phases. 

One way of characterizing the nontrivial properties of 
these phases is the following definition: a state is "triv- 
ial" if it is a product state, or if it can be transformed 
into a state which approximates a product state using a 
local or approximately local quantum circuit. A state is 
nontrivial, or topologically ordered, otherwise. Despite 
the usefulness of this definition in terms of circuits, we 
do not have a similarly useful definition for states at tem- 
perature T > 0. 

Here, we propose an analogous "circuit definition" for 
topological order at T > 0. Another method used pre- 
viously to study topological order at T > is topologi- 
cal mutual information^, a generalization of topological 
entanglement entropy^. While a nonzero topological 
mutual information for T near zero is a useful numeri- 
cal signature of nontrivial phases^, the use of topological 
entropy can give different answers from the circuit defini- 
tion even for at T = 0i£, so both definitions are useful to 
fully characterize a state. Also, it has been shown that in 
many two-dimensional theories there exist stringlike log- 
ical operator s 9 ' 10 . If these logical operators can be gen- 
erated by dragging defects with energy 0(1), then, since 
there is a nonzero density of such defects at T > 0, the in- 
formation which would be topologically protected at zero 
temperature will decohere rapidly at T > (though we 
emphasize that the defect density is exponentially small 
in 1/T; as an experimental example, fractional Hall con- 
ductance is accurately quantized^ even at T > 0) . While 
this kind of "operational definition" is practically useful, 
the circuit definition here will lead to similar operational 
results for all Hamiltonians which are a sum of local com- 
muting terms. 

Conversely, we will show that the four dimensional 
toric codei^ has topological order under the circuit def- 
inition for sufficiently small T > (this will not be a 
rigourous proof since we will hcuristically argue for the 



existence of certain operators using previous results, but 
we then prove that the existence of these operators is 
inconsistent with a trivial state). One advantage of a 
definition using the density matrix rather than the exci- 
tation above the ground state is seen in a toy system in 
the Appendix which is topologically trivial at T = but 
non-trivial at some small T > 0. 

Topological Order at Zero Temperature — We start 
with various previous definitions of topological order at 
T = 0. Topologically ordered Hamiltonians according to 
one definition have a ground state degeneracy that de- 
pends upon the topology of the manifold, and the ground 
state subspace obeys a property called the "disk axiom" 
or "TQO-l".MJ4-i6 ) we qua ntify by: 

Definition 1. Let L denote the system size and let P 
denote the projector onto the ground state subspace. Let 
L* denote some length smaller than L. Then, the ground 
state subspace has (L*,e) topological degeneracy^ if for 
any operator O supported on a set of diameter smaller 
than L* there is a scalar z such that 



\POP- zP\\ < e. 



(1) 



Using circuits we can also define topological order for 
systems with a unique ground state. Consider a unitary 
quantum circuit U where the depth of the circuit multi- 
plied by the maximum range of each unitary in the circuit 
is bounded by some range R. Note that for any operator 
O supported on a set Z, U'OU is supported on the set 
of sites within distance R of Z. 

Definition 2. Let i/jq be the ground state of the Hamil- 
tonian H . We say that the state ipo *s (R, c) trivial if 
there exists a unitary quantum circuit U with range R 
such that \ipQ — Uipprodl < e for some product state tp pr od- 

Certainly, every state is (L, 0) trivial, so we are only 
interested in the case of R < L. Colloquiually speak- 
ing a system will be trivial if it is (R, e) trivial for some 
e << 1 and for some R « L. We can relate the two 
definitions^: if a system has (L* , e) topological degen- 
eracy, then no state ipo hi the ground state subspace is 
(R, S) trivial for R < L*/2 and sufficiently small S, e. To 
see this, suppose ipo is (R, S) trivial. Then, the expecta- 
tion of any operator O in state ip pro d is close to the ex- 
pectation of (ifjo\UOW\ipo). If O is supported on a single 
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site, UOW is supported on a set of diameter less than 
L* and so for such O, (ip pro d\0\ijj pro d) « (U'ViMV'Vi)- 
That is, the subspace spanned by il> pro d and U'ipi has 
(1, 77) topological degeneracy for 77 = 0(e + S). However, 
for 77 sufficiently small compared to inverse system size, 
no such state U^ipi exists since ip pr od is a product state. 

Topological Order at Nonzero Temperature — We be- 
gin by defining a "classical state of range i?" p c i (this 
will replace the use of a product state when we de- 
fine T > topological order) to be a state such that 
p c i = Z~ l cxp(— H c i) where Z is a normalization factor 
and H r j is a Hamiltonian which is a sum of terms all act- 
ing on sets of diameter at most R and all of which are 
diagonal in a product basis. 

Definition 3. A density matrix p is (R, e) trivial if it is 
possible to tensor in additional degrees of freedom ICi on 
each site, defining an enlarged space with Hilbert space 
Hi <8> /Q on each site, such that 

\p-Tr {ICi} (Up cl U^\<e, (2) 

where the \...\ denotes the trace norm (the trace norm of 
a Hermitian operator is the sum of the absolute values 
of its eigenvalues ), where U is a unitary quantum circuit 
with range R, and p c i is a classical state of range R (both 
U and pd are defined on the enlarged space), and where 
the trace is over the added degrees of freedom ICi. 

If we allow H c i to be unbounded and we allow the 
dimension of ICi to be unbounded, then this definition 
is equivalent to saying that p is (i?, e) trivial if it is, up 
to error e in trace norm, equal to an incoherent sum of 
(R, 0) trivial states: p = Y,a P( y a )\^i"v{a)){tj} trlv {a)\, for 
some probability distribution P(a). However, for some 
purposes one might want to construct p using a bounded 
H c i and a bounded dimension on ICi. 

Absence of Topological Order For Two Dimensional 
Hamiltonians With Commuting Terms — We now show 
absence of topological order for any T > 0, under the 
above definition, for two dimensional Hamiltonians which 
are a sum of commuting projectors. The proof is based 
on showing that the density matrix can be approximately 
written as a weighted sum over density matrices of a sys- 
tem with "holes" in it as explained below, and then using 
results fromi^ to write each such density matrix as a triv- 
ial state. 

Consider a two dimensional Hamiltonian H = 
Y2x Qx.i where the terms Qx are commuting projectors. 
Assume that the terms in the Hamiltonian are local, so 
that each projector Qx is supported on some set X which 
has diameter Ri nt which is 0(1). Further suppose that 
each site is in at most 0(1) of such sets X, and for sim- 
plicity consider a square lattice. 

The density matrix is p = Z~ x exp(—/3H), where Z is a 
normalization and j3 = T _1 . Note that for any projector 

Qx wehaveexp(-^Qx) = S Sx£ {o,i}(( 1 " s ^) ex P("ft- 

/ + sx(l — ex p(~/3)) • {I — Qx) ) where we introduce an 
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FIG. 1: Illustration of holes in the lattice. The solid lines indi- 
cate the division of the lattice into squares of size lp, and only 
some of the squares are illustrated. Solid circles indicate holes 
in the lattice, with one hole per square. Dashed lines connect 
holes (dashed line extending outside the solid lines connecting 
to only one hole are intended to indicate connections to holes 
in other square which are not illustrated). Symbols 0, 1, 2, ... 
indicate different subsets of the lattice with holes in it, with 
the dashed lines indicating divisions between those subsets. 

additional variable sx, and sum over sx = 0, 1 (this vari- 
able sx is unrelated to any local spin degrees of freedom 
of the Hamiltonian). Thus, 

p= J2^)lp({s x })P({sx}), (3) 

M 

where the sum is over a set of variables sx, each vari- 
able taking values or 1, with P{{sx}) = — 
Sx) exp(— /?) + sx{l — cxp(— /?)), and where 

p({sx}) = Zdsxjy'H^l-sxV+sxil-Qx)), (4) 
x 

with Z({,sx}) _1 being a normalization. Note that 
p({sx}) is maximally mixed on the ground state sub- 
space of H({s x }) = J2x s xQx- 

We will show that the density matrix p, for any given /3, 
is dominated by a sum over choices of sx in which there 
are lots of "holes" in the lattice, where a hole corresponds 
to a disk Y with radius greater than Ri„t, such that for 
any X with XnY ^ 0, we have s x = 0. See Fig. 1. (Note 
that the variables sx for X not intersecting such a disk 
Y arc not determined by the choice of holes, these vari- 
ables may be either or 1). We divide the square lattice 
into large squares of linear size lp, with lp exponentially 
large in (3 as given in Eq. ([SJl and logarithmically large 
in system size. We call a configuration "valid" if there 
is at least one hole per square. We will show that for 
valid configurations the density matrix p({sx}) can be 
expressed as a local unitary with range R = 21 p acting 
on a classical state and we bound the contribution to p 
from invalid configurations. Combining these results im- 
plies that the sum over sx can be approximated by a 
local unitary acting on a classical state (the error e arises 
from the invalid configurations of sx). 

Consider a given valid choice of sx- For such a valid 
choice of sx, we draw an image such as in Fig. 1. As 
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mentioned, the terms sx may be zero even if X docs 
not intersect a hole. If a given choice of Sx has more 
than one hole per square, we only indicate one hole in 
the square, choosing which one to take according to any 
arbitrary rule (for example, the hole closest to the top 
left corner). We draw dashed lines connecting the holes. 
The dashed lines break the lattice outside the holes into 
regions labelled 0,1,2,... as shown. Then, H({sx}) can 
be re-written as 

H({sx})= J2 H a , b + Y,H a , (5) 

<a,b> a 

where a and b label regions bounded by dashed lines, the 
sum is over neighboring a, b, and R a ,b is supported on 
regions a and b (define H ao = -ffb, a ) and H a is supported 
on region a. 

So, by coarse-graining the lattice, the Hamiltonian is 
a sum of commuting terms, each acting on at most two 
regions. For any a, the terms -ff a ,b and H a which are sup- 
ported on that region all commute with each other. This 
allows us to find a decomposition of the Hilbert space 
which greatly simplifies the description of the problem. 
Eqs. ()6I7[) below differ slightly from lemma 8 o&, but can 
be proven using the same tools of interaction algebra^ 
(see also 2 ^): decompose H a ^ as a sum of product op- 
erators H ao = J2-y^a b (l)Ob b (l)i where the operators 
0° & (7), O b b (j) are supported on a, b respectively and the 
operators 0% (7) are chosen from an orthonormal basis. 
Then, [Of (S), 0^(7)] = for b ^ c, for all 8, 7. Let A ab 
be the algebra generated by the set of 0^(7) for given 
b. The algebras A ab ,A ac commute for b ^ c. 

Let Ha denote the Hilbert space on region a. One way 
for the two algebras to commute is simply that T-L a de- 
composes into a tensor product of Hilbert spaces, and 
each A ah acts on a different space. However, this is not 
the only possibility. Suppose, for example, the Hamil- 
tonians H a ,b are all diagonal in some product basis. For 
example, consider an Ising Hamiltonian with all terms in- 
volving only operators S z . Then all the R a .b would com- 
mute, but we would not have this tensor product decom- 
position. However, we can decompose T~L a into a direct 
sum of Hilbert spaces T-La^ a \ and then further decompose 
each such Hilbert space into a tensor product of 

n ,a(a) 

spaces n a _^ b giving 

n a = 0K (Q) = 0(Ki a) ® (g) K$), (6) 

a(a) o-(o-) <b,a> 

where the product is over b that neighbor a, such that 
each operator H a j, can be decomposed as 

Ha*= £ F fl aW Pf (6) <; W(6) , (7) 

a(a),fl(b) 

where Pa is the operator on H a which projects onto 
UZ [a) and Hf b a) ' m acts on the subspace oiU^ a) 

given by K^®^- 



Define Qa.b to project onto the nonzero energy states 
of H a b, so that the maximally mixed state on the ground 
state subspace of H({sx} is Z^sx}) -1 Yl <a b> (I — 
Qa.b) lla(^ — Qa)i where Z({sx}) is a normalization fac- 
tor. This state is a (2,0) trivial state on the coarse- 
grained lattice of regions because for any choice of the 
variables a(a) for each region a the projection of p({sx}) 
onto the product of spaces Ha is a product state on 
the spaces H"^l'S>'H b ^ a . Thus it is an (i?, 0) trivial state 
on the original lattice for R = 2lp. 

We now bound the contribution of invalid configura- 
tions. Note that given any two sequences {sx } and {s' x } 
such that s' x > sx for all X we have Z({s' x }) < Z({sx}- 
Imagine breaking each large square of linear size lp into 
small squares of size 2Ri nt on each side (this will not 
give the best estimates but simplifies the proof), and 
let nc be the maximum number of projectors that in- 
tersect any of those small squares. Then, the sum of 
Z({sx})P({sx}) over all {sx} such that there is a given 
configuration of holes in those small squares with a total 
of k such holes is at least equal to exp(—(2Rf nt kf3) times 
the sum of the same quantity over configurations with no 
holes. So, the sum of Z~ 1 Z({sx })P({sx}) over configu- 
rations with no holes in any given square is bounded by 

(1 - exp(-R 2 in tP)) (l < i/Rmt)2 \ as follows by counting the 
number of configurations with holes and without. 
Pick 

lp = exp((2i? mt ) 2 /?)i? mt \og{V)/ log(e), (8) 

so (1 -exp(-(2 J R mt ) 2 /3)(^/ i? -*) 2 < e/V. Then, the con- 
tribution of configurations such that at least one large 
square has no holes to the trace in Eq. ([3]) is at most e. 
Thus, for such an lp we can ignore such configurations 
and restrict to a sum over valid configurations giving an 
approximation to p with error at most e in trace norm, 
describing p as an incoherent sum over trivial states; if 
desired, one can write this sum in the form of Eq. ([2]). 

Operational Properties — Consider first T = 0. The 
four dimensional toric code on a torus has a degenerate 
ground state. There are surface operators U x and U z 
that act like the Pauli operators a x and a z on a two- 
dimensional subspace of the ground state space. These 
operators anti-commute so that {U X ,U Z } = 0. Further, 
defining U' x and U' z to denote the same operators trans- 
lated a distance L/2 perpendicular to the given surfaces, 
the ground state expectation value of U'JJ], equals 1, and 
similarly for U' Z U\. At nonzero temperature, the possibil- 
ity of correcting errors^ implies that we can "thicken"— 
those surfaces, giving unitary operators V x and V z which 
are supported within some distance (say, L/8) of the 
given surface which have similar properties at T > 
to the operator U x and U z at T = 0. In particular, we 
expect that for sufficiently small T 

tv{ P V x V z vM) ~ -1) (9) 
tr( P yX f ) » tr(^t) « 1. (10) 
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However, Eqs. (|9I10[) are inconsistent with having an 
(R, e) trivial state for R sufficiently small compared to L 
and e of order unity: if p is an incoherent mixture of states 
iP triv (a), tr(pV^) = Y. a P{o)^tnv{a)KV^ triv {a)). 
Then, since the separation between the supports 
of V x and Vj. is greater than R, this equals 

J2a P ( a ) (^triv (a) | Vj. | Iptriv (a) } (iptriv (a) I V^\lp t riv («)}• 

For this sum to be close to 1 as in Eq. (JTUJ), a iptriv 
chosen at random must be, with probability close to 1, 
an approximate eigenstate of both V x and V z . However, 
this is inconsistent with Eq. ©. One can derive a 
similar inconsistency result for certain defect creation 
processes similar to those in22. 

Conversely, we can comment on this usefulness of a 
trivial state for storing quantum information. To manip- 
ulate quantum information at T > 0, we need operators 
V x and V y which act like the Pauli matrices on a single 
qubit, as in Eq. ©. However, since Eqs. (|9I10|) are incon- 
sistent with having a trivial state, we cannot have also 
corresponding operators V,j,,Vy. Further, such a trivial 
mixed state can be no more useful for storing informa- 
tion than a trivial pure state. Note that it seems too 
much to hope for the converse statement (that a nontriv- 
ial state is useful for storing quantum information), since 
even at T = the Chern insulator— provides an exam- 



ple of a circuit nontrivial state that has a unique ground 
state and so cannot be a quantum memory. 

Discussion — We have proposed a definition of topolog- 
ical order T > 0. While this definition is simple, it allows 
us to make precise statements about how quantum infor- 
mation can be manipulated in two dimensional Hamilto- 
nians which are sums of commuting terms. This raises 
the question of what can happen for arbitrary two dimen- 
sional Hamiltonians. In three dimensions, we expect that 
discrete gauge theories are topologically trivial at T > 
under the circuit definition (see Appendix) . A much more 
interesting question is whether Haah's cod o 24 i 25 , which 
avoids this gauge theory paradigm, is trivial or not at 
T > 0. Using quasi-adiabatic continuation 2 ^, we can re- 
late the circuit definition at T = to whether or not 
one can deform one local Hamiltonian into another while 
avoiding a phase transition; however, an analogous ques- 
tion for T > is open (whether the absence of a phase 
transition in thermodynamic quantities when lowering 
the temperature from infinity to some finite 7/ implies 
that the density matrix at Tf is trivial). Finally, we ask 
if there is a circuit definition for exotic critical points^. 

Acknowledgments — I thank D. Poulin for useful com- 
ments. 
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Appendix A: An Example Hamiltonian With 
Topological Order Only at T > 

Consider the following Hamiltonian. It is defined on a 
four dimensional lattice. On each plaquettc of the lattice, 
we define a spin- 1/2 degree of freedom. On each site, we 
define a degree of freedom with 3 different states, labelled 
0),|1),|2). Define an interaction term for the site degrees 
of freedom, H s , as follows: 

H s = -JJ2 |0M0|®|0>,<0| (Al) 

<i,j> 

-J J2 + ®|2)<(2|) ® (|1},(1| + ®|2) 3 -(2|) 

<i,j> 

+h y £\o) i (o\. 

i 

This Hamiltonian is classical, in that it is diagonal in 
the given basis. We take J > 0. The sums are over 
nearest neighbor The first term favors having on 
site i neighboring on site j. The second term favors 
having 1 or 2 on site i neighboring 1 or 2 on site j. The 
third term favors having on each site for h < and 
favors 1 or 2 for h > 0. The ground state has all sites 
in the state, but one can tune the interaction terms to 
produce a transition at non-zero temperature to a state 
where most spins are in the 1 or 2 state (since one can 
choose each spin to be in either the 1 or 2 state, there 
is an increase in entropy which can outweigh the h term 
in the Hamiltonian), with just a few excitations to the 
state (an isolated is penalized in this state, since most 
of its neighbors arc 1 or 2; it is only when many are 
nearby that one gains energy). The transition is likely 
first order. 

Now add interaction terms for the plaquette degrees 
of freedom. The toric code in four dimensions has two 
terms, which we call A and B, which act on plaqucttcs 
cither surrounding a link or in a cube. These terms are 
diagonal in either the X or Z basis respectively. Define 
new interactions terms C, D as follows: C is equal to A 
multiplied by the product, over sites on that link, of the 
projectors onto the space spanned by the 1 and 2 states 
on those sites and D is equal to B multiplied by the 
product, over sites in that cube, of the projectors onto 
the space spanned by the 1 and 2 states on those sites. 
Finally, we add an additional term E which is equal to the 
sum over plaqucttcs of a z on that plaquettc multiplied 
by the product, over sites bordering that plaquette, of 
the projector onto the the state on those sites. 

Then, we define the Hamiltonian 



H s + C + D + E. 



(A2) 



By tuning h, we can achieve the following phase diagram 
as a function of temperature: the ground state has 
states on all sites. Then, the C, D terms are equal to 
zero, and the E term polarizes the plaquette degrees of 
freedom in the Z direction. Hence, the ground state is a 
trivial product state. 



Trivial 



Trivial 

Topological 



FIG. 2: Illustration of phase diagram expected for the toy 
model. The topological phase is the topologically ordered 
phase of the four dimensional toric code. The trivial phases 
are different trivial phases. The solid line indicates a first 
order transition. The dashed line indicates a phase transition 
in which the site degrees of freedom remain primarily in the 
|1) or 1 2) states, but the toric code portion of the system 
becomes topologically trivial. Scales of axes are arbitrary, 
as is the particular value of h at T = at which the phase 
transition occurs. The point of the model system is that a 
point at T = with h slightly smaller than the critical h can 
change from trivial to topological as temperature is increased. 



However, as we increase the temperature, this changes. 
Consider the state which is maximally mixed over 1 and 
2 degrees of freedom on the sites, multiplied by the toric 
code ground state. This state gains entropy (due to the 
mixture of 1 or 2 states). This state has a different en- 
ergy from the ground state. Part of the difference arises 
from the differing h term, and part of the difference is 
due to the different energy of the expectation value of 
C + D + E in those two states. However, by tuning h 
we can make this state have very slightly higher energy 
density than the ground state. Thus, it seems very likely 
that we can tune h to produce a first order phase transi- 
tion from the ground state at zero temperature to a state 
at some small T > which contains the four dimensional 
toric code. Now, this state at non-zero temperature is 
not exactly the same as the toric code at non-zero tem- 
perature, since any states on the sites will change the 
interaction terms in the toric code Hamiltonian nearby 
leading to some correlation between defects in the toric 
code state; however since those states are rare at low 
temperature (by tuning h we can accomplish the phase 
transition between states at a very low temperature to 
keep the states rare), it is likely that the effect of these 
states will not destroy the topological order in the toric 
code Hamiltonian. 

Thus, while this system is topologically ordered at such 
temperatures under the circuit definition, any reference 
to ground state excitations cannot capture this. There 
is a range of temperatures where one can store quantum 
information in a robust matter, but it is not possible at 
zero temperature for this h. Error correction is accom- 
plished by using the same error correction procedure as 
in the four dimensional toric code. Of course, this exam- 
ple is quite artificial, but other examples may exist close 
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to zero temperature phase transitions: if one is close to a 
first order phase transition at zero temperature, as we are 
here with the phase transition being driven by changing 
h, then increasing temperature can lead to a situation 
where the situation at a nonzero temperature for a given 
value of h (a value that favors the state at zero tem- 
perature) is in the same phase as the zero temperature 
state at a different value of h (a value that favors the 
1 or 2 states). The figure illustrates the phase diagram 
for this model; more realistic models with similar phase 
diagrams probably exist. 



Appendix B: Triviality of Three Dimensional Toric 
Code at T > 

We now show the triviality of the three dimensional 
toric code at T > 0. It is possible that this method could 
be extended to show the triviality of other discrete gauge 
theories in three dimensions; our approach relies on con- 
structing operators with certain commutation relations 
but probably similar operators exist for other discrete 
gauge theories. 

The three dimensional toric code has topological order 
at T = 0. It has a phase transition at non-zero tempera- 
ture, and the topological entropy becomes non-vanishing 
below this phase transition temperature^. 

The Hamiltonian is defined by a three dimensional cu- 
bic lattice, with spin- 1/2 degrees of freedom on the links 
of the lattice. We define two operators, plaqucttc and 
star operators, by 

b p =n°?> ( bi ) 

iEp 

A s =l[af. (B2) 

These operators commute and the Hamiltonian is denned 
by 

H = -\ a Y j A s -\bY j B v ( B3 ) 

s p 

Much of the three dimensional case follows the two 
dimensional case. Similar to Fig. 1 in two dimensions, in 
three dimensions we divide the cubic lattice into cubes 
of linear size lp. We again write the thermal density 
matrix p as a sum over p({sx}), as in Eq. Q. As in two 
dimensions, for £g taken to be of a size which is of order 
log(V)/log(e), one can show that the sum is dominated 
by "valid" configurations, where now a configuration is 
valid if there is at least one hole per cube. 

However, at this point the three dimensional result 
needs new techniques. In two dimensions, we showed 
that each p({sx}) for any valid configuration was trivial 
using techniques based on ideas fromi^ which work for 



any Hamiltonian which is a sum of commuting projec- 
tors. However, in three dimensions, we need to do some- 
thing different, since the Hamiltonian H({sx}) cannot 
be coarse-grained to a Hamiltonian in which each term 
acts on only at most two sites (we would need to have 
long "line-like" holes to be able to do this whie the holes 
are only small "point-like" configurations here). So, we 
develop a technique that works for the three dimensional 
toric code to show that p({sx}) is trivial for valid config- 
urations. This method relies on constructing particular 
operators using Pauli spin operators of the toric code. 

Consider the operator that we call Fj jS , where i labels 
a bond and s labels a site that neighbors the bond i. This 
operator is defined by 

Fi, s = erf A s . (B4) 

That is, Fi tS is the product of ia y on site i with crj over 
sites j i with j G s. We consider, as mentioned, the 
case that s neighbors i. Thus, there are two possible 
choices of s for any given i, and those choices can be 
thought of as defining a direction: s is defined by the 
bond i and the direction that one must move on the bond 
i to arrive at site s. 

Note that F i s commutes with all operators B p . Also, 
Fi_ s commutes with all operators A t except that it anti- 
commutes with the two operators A t with t neighboring i 
(one of these operators A t is equal to A s while the other 
is obtained by moving on the opposite direction on bond 
i). Define the unitary 

Ui, s = exp(i^F itS ). (B5) 

Because of the commutation relations above, this unitary 
obeys: 

U^BpUl = B p , (B6) 

and 

i^t ^ U hS A t ul = A t , (B7) 

where the notation i ^ t means that bond i is not a 
neighbor of site t. However, 

U^AsUl = of. (B8) 
We consider a general class of Hamiltonians 

H(S, V, B) = -X A J2 A, - X b B v ~ *A J2 ^ ( B9 ) 

sG5 p£V i£B 

where S, V, B are sets of vertices, plaquettes, and bonds, 
respectively. Thus, the three dimensional toric code 
Hamiltonian H is H(S,P,B), with S,V being the sets 
of all vertices and plaquettes respectively, while B is the 
empty set. 

Every configuration of the {sx} defines a Hamiltonian 
H(S,V,B), with B being the empty set, and S,V being 
some subset of the vertices and plaquettes, depending 
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upon which sx are non-zero. Now, consider a configu- 
ration of the {sx } and the corresponding Hamiltonian 
H(S,V,B). 

We define a certain unitary transformation for each 
closed surface on the dual lattice. A surface on the 
dual lattice is made of plaquettcs on the dual lattice and 
each such plaquette intersects some spin- 1/2 degree of 
freedom. We assume our surface is orientable (this will 
hold for all the surfaces described below, which are con- 
tractible, but will exclude perhaps certain exotic situ- 
ations on certain three dimensional manifolds), and we 
imagine an arrow on each bond which intersects the sur- 
face; this arrow will point in the direction of the orienta- 
tion. For each such surface, and each bond i, let s(i, £) be 
the vertex reached by starting on bond i and moving in 
the direction of that arrow. We will be interested in oper- 
ators Ui jS (i^-\ below. Let s(i, £) be the vertex reached by 
starting on bond i and moving opposite to the direction 
of the arrow. 

Let us say that this surface is a "free surface" if for all 
bonds i which intersect the surface, the operator A g , t ^ 
does not appear in the Hamiltonian H(S,V,B) and if 
none of the operators F i s ^ ^ intersect any bond in B. 
Note that if sx contains a hole, then the surface sur- 
rounding that hole, oriented outwards, is a free surface 
for H(S 7 V 7 B). 

Take any valid Hamiltonian H({sx}), take any free 
surface £ surrounding a hole and orient that surface 
outwards. Take any operator for any i inter- 

sected by the free surface. Consider the Hamiltonian 
Ui lS (i t £)H({sx})Uj This define some new Hamil- 

tonian; the plaquette terms B p are unchanged by this 



producedure, as all are all the vertex terms A s except for 
the vertex term A s for s = s(i,£) which is turned into the 
operator erf by Eq. (jB8|) . Note that since the surface is 
free, there are no A s intersecting the surface on the inside 
of the hole (such an operator occurs for s = s(i, £) which 
does not appear by definition of a free surface); the only 
such A s are those on the outside of the hole. Note that 
this new Hamiltonian Ui tS (i t g)H ({sx})Uj a ^ ^ has a new 
free surface: this free surface is given by taking £ and 
moving it outwards along bond i. Repeating this proce- 
dure, we can move each free surface outwards. Note that 
since the original Hamiltonian has B = 0, all operators 
erf that are produced by this procedure are produced far 
enough inside the free surfaces that they do not stop the 
surfaces considered in future steps of this procedure from 
being free; that is, such operators are never transformed 
by the operators Ui tS u£\ in later steps. 

We then choose a specific order to apply this proce- 
dure: in the first round, we take each free surface sur- 
rounding each hole, and iterate this procedure until that 
free surfaces reach the boundary between cubes. This 
first round (which is the first round of the quantum cir- 
cuit that we use to show that p({sx}) is trivial) is a 
product of unitaries, one on each cube, each such uni- 
tary being the product of the Uigug) used to move that 
surface out until it hits the boundary. Then, in the next 
round, we move the free surfaces one further step out- 
wards, so that they intersect; once the free surfaces in- 
tersect, the terms A s are removed from the Hamiltonian. 
The final Hamiltonian contains only terms diagonal in 
the S z basis, and hence is classical. Thus, the original 
Hamiltonian H({sx}) is trivial. 



